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Abstract: Investigating the production of polarized heavy quarkonia in terms of recently 
proposed QCD factorization formalism requires the knowledge of a large number of input 
fragmentation functions (FFs) from a single parton or a heavy quark-antiquark pair to a 
polarized heavy quarkonium. We study these universal FFs at the input factorization scale 
pLQ > 2mQ, with heavy quark mass mg, in the framework of nonrelativistic QCD (NRQCD) 
factorization. We express these FFs in terms of perturbatively calculable coefficients for 
producing a heavy quark-antiquark pair in all possible NRQCD states, multiplied by corre¬ 
sponding NRQCD long-distance matrix elements for the pair to transmnte into a polarized 
heavy qnarkoninm. We derive all relevant NRQCD operators for the long-distance matrix 
elements based on symmetries, and introdnce a self-consistent scheme to define them in ar¬ 
bitrary d-dimensions. We compnte, np to the first non-trivial order in Ug, the pertnrbative 
coefficients for prodncing a heavy qnark pair in all possible S-wave and P-wave NRQCD 
states. We also discnss the role of the polarized FFs in generating QCD predictions for the 
polarization oi J/fj prodnced at collider energies. 



Contents 


1 Introduction 2 

2 Heavy quarkonium polarization in NRQCD factorization 5 

2.1 General selection rules 5 

2.2 Nonvanishing channels for the polar angular distribution 7 

2.3 NRQCD power counting 9 

3 Definition of polarized NRQCD LDMEs in d dimensions 11 

4 NRQCD factorization on input FFs 13 

5 Discussion and Summary 15 

A Polarized NRQCD LDMFs 16 

A.l NRQCD LDMEs 16 

A.2 Projection Operators 18 

A. 3 Expand the polarized LDMEs to NLO in perturbative NRQCD 19 

B Single-Parton Fragmentation Functions to a Polarized Heavy Quarko¬ 
nium 20 

B. l Gluon FFs 20 

B.2 Different quark FFs 21 

B.3 Same quark FFs 21 

B.l ^Pj operators with orbital angular momentum summed 22 

B. 5 Comparing with other calculations in the literature 22 

C Double-Parton Fragmentation Functions to a Polarized Heavy Quarko¬ 
nium 23 

C. l A-functions 23 

C.2 Leading Order 24 

C.3 NLO - Vector 24 

C.l NLO - Axial-vector 25 

C.5 NLO - Tensor 26 

C.6 ^Pj operators with orbital angular momentum summed 27 

C.6.1 Leading Order 27 

C.6.2 Next-to-leading order 28 


- 1 - 




1 Introduction 


Since the discovery of the J/il>, heavy quarkonia, with their clearly separated multiple mo¬ 
mentum scales, have been serving as ideal systems to test our understanding of QCD bound 
states and their hadronization processes. Unfortunately, a theoretically and phenomenolog¬ 
ically satisfying framework is still lacking for heavy quarkonium production. The problem 
is more acute with the recently discovered XTZ-mesons, since the production of these ex¬ 
otic mesons requires good understanding of the production of conventional heavy quarkonia 
[ 1 , 2 ]. 

A phenomenologically successful model for heavy quarkonium production is based on 
the non-relativistic QCD (NRQCD) factorization [3, 4], which factorizes the production 
cross section into the production of a heavy quark pair, multiplied by the transition for 
the pair to transmute into the observed quarkonium. The production of the heavy quark 
pair is effectively perturbative and is organized in powers of Ug and v, the relative velocity 
of the heavy quark in the pair’s rest frame; and the corresponding transition is nonper- 
turbative and is represented by a set of NRQCD long-distance matrix elements (LDMEs). 
If the factorization is correct to all orders in these LDMEs will be universal, i.e. pro¬ 
cess independent. For heavy quarkonium production, although a formal proof of NRQCD 
factorization to all orders in Og is still lacking, it has been shown to be valid at the next- 
to-leading order (NLO) in many processes. With certain modification, it also works at the 
next-to-next-to-leading order (NNLO) in some specific cases [5-7]. Phenomenologically, 
NLO NRQCD factorization calculations successfully explain the transverse momentum px 
distribution of many heavy quarkonium states with a few NRQCD LDMEs fitted from the 
data [1, 2]. 

Nevertheless, the current NRQCD factorization formalism is far from perfect in describ¬ 
ing data on the heavy quarkonium production at high pT- It is understood now that the 
fixed-order NRQCD calculation suffers from large high-order corrections, due to the large 
power enhancement in the forms of p'^/niQ and the large \og{p^/m^)-iype logarithms from 
high orders in Ug [8], where ruQ is the mass of heavy quark. For example, the NLO NRQCD 
calculation for the yield of J/'0 is orders of magnitude larger than the leading order (LO) 
calculation for and ^Pf'^ channels at px > 15 GeV [9-12]. Although the most phe¬ 
nomenologically important leading power contribution is claimed to be included in the NLO 
NRQCD calculation [11, 13], the existence of large logarithm may potentially undermine 
the convergence of ag expansion at large pT ^ mg. 

Recently, a new QCD factorization formalism has been proposed to study heavy quarko¬ 
nium production at large pt [8, 14-18]. In this formalism, the cross section is expanded 
by powers of mq/p^. It was proved to all orders in Ug that the dominant leading-power 
(LP) terms, as well as the next-to-leading power (NLP) terms, can be systematically fac¬ 
torized into the perturbatively calculable hard parts for producing a single parton (or a 
heavy quark pair at NLP), convoluted with corresponding single parton (or heavy quark 
pair) fragmentation functions (EEs) to the observed heavy quarkonium. All nonpertur- 
bative contributions are included in these EEs, whose scale dependence is determined by 
a closed set of evolution equations with perturbatively calculable evolution kernels. By 
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solving the evolution equations, large perturbative log(p|./mQ)-type logarithms could be 
resummed to all orders in Og. Because of the systematic treatment of the powers of P^/ttIq 
and log(p^/mQ), the QCD factorization is expected to converge faster in Ug expansion than 
the NRQCD factorization. 

With the LO partonic hard parts for hadronic production of heavy quarkonia at high px 
calculated in Ref. [18], and the evolution kernels of FFs available in Refs. [8, 16, 17], a set of 
FFs at some input scale po > 2mQ is needed and necessary for evaluating the rate of heavy 
quarkonium production, and for any phenomenological application of the QCD factorization 
formalism. These input FFs are only unknowns in the QCD factorization framework, and 
are sensitive to the observed heavy quarkonium states and their polarizations. In principle, 
these input FFs should be extracted from data, like the case for pion and kaon production at 
high pt- Nonetheless, this task is extremely hard in practice, especially with the inclusion of 
NLP contributions. For example, it requires at least four one-variable single parton FFs and 
six three-variable heavy quark pair FFs for evaluating the production rate of polarization- 
summed J/il] [19, 20]. The number of FFs should double for calculating the cross section 
of producing polarized J/'i/’- 

However, different from the pion and kaon FFs, the heavy quarkonium FFs have an 
intrinsic large scale mg 3> AqcD) thus they are partially perturbative. As early as in 1993, 
Braaten, Cheung and Yuan proposed to apply the Color-Singlet Model and the NRQCD 
factorization to further separate the perturbative and nonperturbative contributions to the 
input FFs [21, 22]. By choosing the input QCD factorization scale pQ > 2mQ and NRQCD 
factorization scale p\ ~ mg, neither PQ/rnq nor mg/p\ is large, and consequently, the 
NRQCD expansion of the input FFs is expected to have a fast convergence. In Refs. [19, 
20], we calculated the polarization-summed input FFs in the framework of the NRQCD 
factorization. We derived, up to the NLO in as, the perturbative hard parts for an energetic 
single parton (or a heavy quark pair) to fragment into a heavy quark pair in all possible 
S-wave and P-wave NRQCD states, while the NRQCD LDMEs cover the nonperturbative 
transition rate for the pair to transmute into the observed heavy quarkonium. ^ These input 
FFs are needed for evaluating the factorized LP and NLP contributions to the production 
of polarization-summed heavy quarkonia. 

With our calculated input FFs, the first test of the QCD factorization formalism for 
heavy quarkonium production is shown in Ref. [23]. It was found that without the re¬ 
summation of large logarithms, the simple and fully analytical LO calculation in the QCD 
factorization approach successfully reproduces the very complicated and numerical NLO 
NRQCD calculations for the yield of J/i!) at px ^ 15 GeV, for all phenomenologically im¬ 
portant and relevant channels. It is very important to note that the NLP contribution, 
although suppressed by a factor of rug/p^ compared to the LP contribution, are crucial for 
and channels even when px approaches 100 GeV. A similar calculation with only 
the LP contribution cannot reproduce these two channels [24]. This finding clearly demon¬ 
strates the importance of the NLP contributions to the heavy quarkonium production. 

^Most single parton input FFs in Refs. [19] are not new. See Refs. [19, 20] for the details and references 
therein. 
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Since the polarization is an important observable for exploring the production mech¬ 
anism of heavy quarkonium [25-27], it is critically important to study the production of 
polarized heavy quarkonia in the framework of the QCD factorization formalism, which 
requires a set of polarized input FFs. In this paper, we extend our work on the polarization- 
summed input FFs [19, 20] to the study of the polarized input FFs by using the NRQCD 
factorization approach. Different from the production of unpolarized heavy quarkonia, we 
show that, if the polarization of the produced heavy quarkonium is observed, many more 
NRQCD channels at the same order in v expansion are needed. Especially, off-diagonal 
channels are needed for the production of some heavy quarkonia even at the leading power 
in V. We demonstrate that conservation laws, angular momentum addition rules, and veloc¬ 
ity scaling rules, as well as the heavy quark symmetry, are very important for constraining 
the size of contributions from various channels. 


Similar to the polarization-summed input FFs [19, 20], the intermediate steps in de¬ 
riving the polarized input FFs require the perturbative treatment of ultraviolet (UV) and 
infrared (IR) divergences. Conventional dimensional regularization (CDR) is the most con¬ 
venient tool to regularize all these divergences. However the application of CDR on the po¬ 
larized input FFs is not straightforward, since the definitions of polarized NRQCD LDMEs 
in an arbitrary d-dimension are very non-trivial. In principle, one needs to consider the 
coupling of two d-dimensional heavy quark spinors, and L-S coupling if the orbital angular 
momentum is nonzero, e.g., for the channels. In this paper, by requiring the symme¬ 
tries to be preserved when generalizing the polarized NRQCD LDMEs from 4-dimensions 
to d-dimensions, we provide a simple scheme to separate the contributions with different J 
and \Jz\. 


The rest of this paper is organized as follows. In Section 2, in terms of the NRQCD fac¬ 
torization formalism, we derive, by using conservation laws and angular momentum addition 
rules, all nonvanishing NRQCD LDMEs for the production of polarized heavy quarkonia. 
We single out the derivation of the nonvanishing NRQCD LDMEs that contribute to polar 
angular distribution of heavy quarkonium production, which are of the most phenomenolog¬ 
ical interest. We also discuss how velocity scaling rules help simplify problems in practice. 
Then in Section 3, we introduce a scheme to dehne these polarized NRQCD LDMEs in an 
arbitrary d-dimension. In Section 4, we outline the technical steps needed for calculating 
the single-parton and heavy quark-antiquark pair {QQ pair) FFs within the NRQCD factor¬ 
ization approach. Since the calculation is very similar to that for the polarization-summed 
case, which has already been explained in great details in Refs. [19, 20], we do not repeat it 
in this paper, but providing all mathematical tools specific to the calculation of polarized 
input FFs in Appendix A. Finally in Section 5, we briefly discuss the impact of our results 
for the polarized J/'0 production. All of our results on the polarized input FFs are listed 
in Appendixes B and C. In Appendix B.5, we compare our full results of single parton FFs 
with some input FFs available in the literature. 
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2 Heavy quarkonium polarization in NRQCD factorization 

In this section, we start from experimental observables to find all allowed NRQCD channels 
for polarized heavy quarkonium production. In Subsection 2.1, we drive the most general 
selection rules for nonvanishing NRQCD channels, by employing only the conservation laws 
and angular momentum addition rules. These selection rules are valid for any observables 
related to the heavy quarkonium polarization. Then in Subsection 2.2, we focus on the 
polar angular distribution of decay products of produced heavy quarkonia, and derive more 
selection rules for this specific observable. Finally in Subsection 2.3, we adopt the NRQCD 
power counting rules to find the relative importance of the survived channels for the polar 
angular distribution. 


2.1 General selection rules 


Experimentally, the polarization of a produced heavy quarkonium is determined by mea¬ 
suring the angular distribution of its decay products in its rest frame. Taking J/iIj as an 
example, the angular distribution of the decaying /+/“ lepton pair can be parameterized as 
(see Ref. [28] for a thorough discussion, and [29] for a recent review) 


- oc 1 -|- Ag cos 9 + sin 9 cos 2^p + sin 29 cos ip 




( 2 . 1 ) 


-|- sin 9 sin 2ip + A^^ sin 29 sin ip, 


where 9 and ip are the polar angle and the azimuthal angle of respect to a chosen z-axis 
of the pair’s rest frame, respectively. 

Similarly, one can do the decomposition for any heavy quarkonium state, such as XcJ 
(see [30-32]). In the narrow-width approximation for producing a heavy quarkonium, all 
the coefficients A’s in Eq. (2.1) and corresponding coefficients for other heavy quarkonium 
states can be related to the heavy quarkonium’s spin density matrix pjjj jjj , where 

{J^) is the z-component of the spin of the heavy quarkonium H in the amplitude (complex 
conjugate of the amplitude). Here the quantum interference with ^ is allowed. 
Comparing theoretical predictions of these A’s with the data serves as an important test of 
the theory [33, 34]. 

Within the framework of the NRQCD factorization [4], the heavy quarkonium’s spin 
density matrix can be factorized as 




E 


= ? 


[QQ{n,n)\ 


J(0|O 


| 0 ), 


( 2 . 2 ) 


where the scale dependence is suppressed, is the spin density matrix of the polarized 

intermediate QQ pair with quantum numbers n = j and h = 

amplitude and the complex conjugate of the amplitude, respectively, where S{S), L{L), 
J(J), and Jz{Jz) are spin, orbital angular momentum, total angular momentum and total 
angular momentum along the z-axis of the pair in the amplitude (complex conjugate of the 
amplitude), respectively. Similarly, 5, 6 = 1 or 8 represent the color-singlet or color-octet 
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state of the produced pair. In Eq. (2.2), the NRQCD LDMEs ’-h^|0) represent 

nonperturbative probabilities for the intermediate [QQ{n,n)] states to hadronize into the 
polarized heavy quarkonium H{J^, J^). In principle, for each channel with the dehnite n 
and n in Eq. (2.2), there are still many LDMEs from contributions at different powers of 
v‘^. In this section, we do not try to explicitly separate different LDMEs contributing to the 
same channel in Eq. (2.2), since we only focus on the selection rules for different channels. 



Figure 1. Schematic picture of angular momentum coupling in the hadronization from QQ-pair 
to heavy quarkonium H. 


The summation in Eq. (2.2) runs over all n and h, but, the LDMEs \Qi) 

vanish for some combinations of n and h. A particular channel [QQ{n, h)] can contribute to 
the spin density matrix only when n and h have the same conserved quantum numbers. For 
example, 6 = 6 and ( — 1)^+^ = (—since the color and parity are always conserved in 
QED and QCD. With color 6 = 6, the state of the heavy quark pair, [QQ[n^n)\ dehned in 
NRQCD, should be an eigenstate of charge conjugation, which leads to (—= (—l)^"’"'^. 
Combining this with the conservation of color and parity, we obtain the following selection 
rules: 


6 = 6, S = S, |L-Z| =0,2,4--- , (2.3) 

for the LDMEs V) m Eq. (2.2). 

Furthermore, angular momentum addition rules can provide more constraints on the 
choice of n and h for the LDMEs ’-L^|0). These nonperturbative NRQCD LDMEs 

[Q/C^(7T.,Tl jj 

can be represented schematically by the diagram in Fig. 1, which is shown in cut diagram 
notation, in which the amplitude and its complex conjugate are combined into a forward 
hadronization diagram for a heavy quark pair to evolve into the heavy quarkonium H and 
unobserved particles X, and the hnal state is identihed by a vertical line. In Fig. 1, is 
the total angular momentum of the heavy quarkonium L6, {J^) is its z-component in 

the amplitude (complex conjugate of the amplitude), and and are the total angular 
momentum and its z-component of the unobserved particles A, respectively, which also 
include the relative orbital angular momentum between H and X. By summing over all 
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possible states of X, the NRQCD LDMEs, as shown in Fig. 1, can be expressed as 


{0\O 


[QQ{n,n)] 


|0) E J, J){J, Jz\A J^){A J^\J, Jz), 

X JXJX 


(2.4) 


where the coefficient A{J^, J, J) depends on the dynamics, but, is independent of the 

z-component of any angular momentum due to the rotational invariance of QCD. Eq. (2.4) 
leads to the following constraints, 

IJ — J| < 2J^, and Jz — Jz = ■ (2.5) 


Eqs. (2.3) and (2.5) summarize all the constraints for the NRQCD LDMEs in Eq. (2.2) for 
the production of a polarized H that decays into two bodies, with the measured full angular 
distribution of the decaying products, like the J/'0-production in Eq. (2.1). 

The constraints in Eq. (2.5) can be modified if the angular distribution of the decay 
products is integrated. As we will show in the next subsection, the polar angle 0-distribution 
(after the integration of the azimuthal angle of the decaying products) depends only on 
the diagonal entries of the heavy quarkonium spin density matrix pjjj Thus, only the 

diagonal LDMEs with = J^, are needed for the 0-distribution of 

heavy quarkonium production, and therefore, the constraints in Eq. (2.5) can be simplified 
to 


J — J| < 2J^, and Jz = Jz-, (for the polar 0-distribution). (2-6) 


For the yield of a polarization-summed heavy quarkonium Lf, we integrate over both the 
polar 0 and azimuthal ip distribution of the decay products, and consequently, we only need 

By choosing and summing over jf in 

Eq. (2.4), the orthonormality condition of angular momentum requires 

J = J, and Jz = Jz, (for producing a polarization-summed H). (2.7) 


That is, for the yield of a polarization-summed H, only the fully diagonal entries in both 
J and Jz of the heavy quarkonium spin density matrix are needed. 


2.2 Nonvanishing channels for the polar angular distribution 


In this paper, we focus on the polar angular distribution, which is of the most phenomeno¬ 
logical interest. For J/il) production, the polar angular distribution means the \q term in 
Eq. (2.1). By integrating the azimuthal angle ip from 0 to 27r on both sides of Eq. (2.1), we 
obtain 




oc 1 -|- Ae cos^ 0. 


( 2 . 8 ) 


dcos 0 

In the narrow-width approximation, the polar angular distribution for producing the J/'0 
that decays into a lepton l~^l~ pair can be expressed as 




E 


jJ/4' 

•J Z z 


J/1/> 

P jJ/i^ fJ/’/' 




(2.9) 
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dcos 9 




where decay is the decay matrix element of J/V’j and €i{Jz^'^) and e|(j/^^) are the polar¬ 
ization vectors of the J/'0 in the decay amplitude and its complex conjugate, respectively. 
Due to the parity conservation and the rotational invariance about the z-axis after the 
integration of azimuthal angle, (p, the decay matrix for the polar angular distribution must 
have the structure ~ + *^2 5^’^6^’^, where ci and C 2 are the coefficients deter¬ 


mined by the dynamics. That is, only the diagonal entries, i.e. the terms with 

of the spin density matrix, in Eq. (2.2), can contribute to the polar angular 

distribution in Eq. (2.9) [29]. 

Furthermore, because of the parity and time-reversal invariance of QED, we have 

Consequently, the polar angular 
distribution depends only on the combination of j/^ + j/^- Similar argu- 

ment also applies to the production of other heavy quarkonium states. That is, the polar 
angular distribution of the decaying products from any produced heavy quarkonium H 
depends only on the combination of the spin density matrices: pjlj ju -t- P^jh -jh- 

Therefore, we can use the selection rules in Eqs. (2.3) and (2.6) to find out all relevant 
QQ-paii channels (n, fi) that contribute to the polar angular distribution. In addition to 
all diagonal channels with n = h, there are also many offdiagonal channels. For example, if 
n is S'-wave, all S-D mixing channels are allowed; and if n is P-wave, there are P-P mixing 
channels as well as P-F mixing channels: 


*^0 > ^ 2,0 
Ipi^] 


l,Jz' 3,Jz 


3c.[6] SnW 


ipW 3p 
^J,Jz^ ^ 


J,Jz 




( 2 . 10 ) 


where 6=1,8 labels the color of the QQ pair. 

Since the polar angular distribution is proportional to the sum of two spin density 
matrices with the opposite signed z-component of quarkonium angular momentum, we 
introduce the reduced spin density matrices as 

Pfj?\ = 2 ^ P-J¥,-J¥) = P^ 

with X = L,T, TT, • • • , corresponding to \J^\ = 0,1, 2, • • •, respectively. From the NRQCD 
factorization in Eq. (2.2), and using V) = ( 0 |C>^twhere -n = 

and similarly for —h, which is an immediate consequence of the time-reversal 
invariance of QCD, we hnd that the reduced spin density matrices in Eq. (2.11) can be 
factorized as. 


P? = E ^ (pLf + PS-n) t ((0|O«P,,„_,„|0) + {0|O"-g,_. _„,||0)) , 


where the reduced NRQCD operators are dehned as 
qHx = ^ 1 

I A" I 


^ ,-jf^ 


( 2 . 12 ) 


where the subscript A is defined right below Eq. (2.11). 
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(2.13) 


Table 1. Essential channels for various polarized heavy quarkonium production, with relative 
power-counting of each channel explicitly. A represents the polarization of the heavy quarkonium, 
and could be L, T or TT depending on the specific state. Q can be either charm or bottom. 


Quarkonium 

Essential channels 

vginS) 

(‘4‘i)b"i. (’d*!)!-’! 

^lj{nS),T{nS) 

(3sSl)b»|. csfi)!..’], (’si'bbi, 

hginP) 

(iFd)[«»i. {'4*')Ki. {=4®1)[<’“i. (=‘s1t)i»‘’i 

xgj{nP) 



2.3 NRQCD power counting 

The general discussions in last subsection were based only on the symmetry properties 
of QED and QCD. In principle, one needs all allowed QQ channels for any quarkonium 
production. In practice, however, different channels are not equally important. The relative 
importance of different channels in power of v is governed by the velocity scaling rule of 
NRQCD [4]. 

To derive the scaling rule, one first expands the wave function of a physical heavy 
quarkonium H into different Fock states in NRQCD effective field theory. For a conventional 
heavy quarkonium, its dominant Fock state is a QQ-pair Fock state with a definite quantum 
number, which is the same as the one usually used to denote the produced quarkonium, 
while all other Fock states are suppressed by powers of v. The relative suppression can be 
estimated by color multipole expansion: the suppression is at 0{v) if two Fock states are 
related by a color electric dipole transition (FI); and the suppression is at if the 

two Fock states are related by a color magnetic dipole transition (Ml) [4]. Therefore, the 
expansion for the wave function of a heavy quarkonium H with quantum number 
can be expressed as, 

\He^+^Lj)) =O(1)|Q0(25+1lW)) + 0{v)\QQe^+\L ± l)[?j)5) ^ 

+ + 0{v^)\QQC^+^L'^i;%g) + • • • . 

In addition to the suppression caused by the Fock states expansion, there is another suppres¬ 
sion factor at for the QQ-pair channel (n,h) because of the derivative operation 

associated with the orbital angular momentum of the operators. Based on these power 
counting rules, one can easily work out the relative importance in v of different channels 
for the production of a given heavy quarkonium. We list the relative power counting of 
phenomenologically important channels for S'-wave and P-wave quarkonium production in 
Table 1. Three general conclusions are in order: 

• First, a channel mixing between different waves is always suppressed by at least 0{v‘^). 

For example, the 

comparing with 


*^0 > ^ 2,0 


0 ’ ^0 


channel in Fq. (2.10) is always suppressed by 0{v‘^) 
channel for any quarkonium state. 
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• Second, a channel mixing between different total angular momentum but the same 
orbital angular momentum can be expressed in terms of a “diagonal” channel up to 
0{v'^) correction. For example, the mixing channels (2-10) 

contributes to the production oi H = J/V’, ip{nS) or T(ns). If we assume the heavy 
quark spin symmetry, the spin of the intermediate QQ-Y>sh in the amplitude, Sz, 
and that in the complex-conjugate amplitude, Sz-, should equal to the spin of the 
hadron . Furthermore, the NRQCD LDMEs corresponding to the mixing channel 
and the channel with Sz = Sz = are equal to each other due to heavy quark 
spin symmetry. As a result, we can dehne a “diagonal” channel ^3>'S’zp[8]^ 3,Szp[8]^ 
study heavy quarkonium production, in which the spin of the QQ-pair is Sz, and the 
orbital angular momentum is summed over. As heavy quark spin symmetry holds up 
to 0{v‘^) correction [4], this approximate treatment is valid at the same order of the 
precision. 


• Third, for any heavy quarkonium state, there is at least one S-wave channel which 
contributes at leading power in v, for any P-wave and higher-wave heavy quarkonium 
state, there is at least one P-wave channel which contributes at leading power in v, 
and so on. Thus lower-wave channels are usually more important. 


With this consideration, we will calculate only S'-wave and P-wave channels in this paper, 
including the “diagonal” channels (^^,Szp[b]^3,Szplb]^ reduced from P-P mixing channels in 
Eq. (2.10), and leave a complete treatment of offdiagonal channels and higher-wave channels 
for future study. Since it involves only diagonal channels in our treatment (up to O(u^) 
corrections), the factorization formula in Eq. (2.12) can be further simplihed to 




; _ ( plQQ] P plQQi ] _ ((01 10) -h (OlO^^- 

n 




(2.15) 


where, similar to the parameter A dehned in Eq. (2.11) for the z-components of heavy 
quarkonium angular momentum , we introduced a parameter, X' = L,T,TT, ■ ■ ■ repre¬ 
senting the z-component of the QQ-paii's angular momentum, \Jz\ = 0,1,2, •••, respec¬ 
tively. For example. 


dQQ] _ 


PriT 2-^ 2 V ^ P-n-n j — ^{n,n) ’ 

n{\JA=^) |Jz|=l 

n{\Jz\=l) 


(2.16) 


| P |=1 


To relate our results in this paper to polarization-summed results, we sum over the polar¬ 
ization A in Eq. (2.15), and obtain 


A 


^E 

np 




EP®' 

ny 


1 


‘A' 


2J + 1 




(2.17) 
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where |0) = is the usual polarization-summed NRQCD 

LDME. To obtain the second equal sign in Eq. (2.17), we have used the factor that the 
summation )]l^^ independent of because of rotation invariance. 

Our calculations in this paper should be sufficient for the current phenomenological 
study of S'-wave and P-wave heavy quarkonium production. In the next section, we provide 
a simple scheme to dehne all these operators in d-dimensions. 


3 Definition of polarized NRQCD LDMEs in d dimensions 

In Refs. [19, 20], we calculated both the single parton and heavy quark pair EEs for an 
polarization-summed heavy quarkonium. We used CDR to regularize all the UV and IR 
divergences. By doing this, we implicitly generalized the polarization-summed NRQCD 
LDMEs to arbitrary d-dimensions. With of the heavy quarkonium summed over, this 
generalization is simple since there is no special direction in the heavy quarkonium rest 
frame. In this polarization-summed case, the possible generalization can be achieved by 
extending the SO{3) symmetry to the SO{d — 1) symmetry [19, 20]. That is, we demand 
all QQ-pair channels to be covariant under the action of SO{d — 1) rotation group, similar 
to the 4-dimensions case. 

The situation is more complicated with the quarkonium’s polarization observed, where 
a specihc z-axis direction needs to be specihed. For example, in the hadron helicity frame, 
the z-axis is chosen to be along the moving direction of the heavy quarkonium in the 
Laboratory frame. To separate contributions with the same J but different jj^j, we need 
to know more details of the SO{d — 1) rotation group. The situation could be more severe 
with angular momentum couplings, such as the L-S coupling, which is exactly what we have 
in the NRQCD factorization. 

In this section, we provide a simple scheme to separate the contributions with different 
J and \Jz\, which is required for the calculation of the polar angle distribution of the decay 
products, like the X$ in Eq. (2.8). In our scheme, we only require: (1) the wave function of 
the heavy quark pair preserves all the symmetries about the z-axis when it is generalized to 
d-dimensions, and (2) the d-dimensional polarization-summed NRQCD LDMEs in Ref. [19, 
20] are recovered after adding up corresponding polarized ones. Consequently, the following 
rules are valid for 5-wave and P-wave channels, 

• Case 1: for ^5i and ^Pi, the wave functions of the heavy quark pair with | J^j = 1 (or 
Jz = 0) in its rest frame are antisymmetric (or symmetric) when flipping the direction 
of all axes except i axe. Note that, it is hard to dehne the rotation operation in d 
dimensions, but the hipping operation described here is always well-dehned^. 

• Case 2: for ^Pj, the wave function with J = 0 is a scalar, i.e. it has a SO{d — 1) 
symmetry. Wave functions with J = 2 (or J = 1) are symmetric (or anti-symmetric) 
tensors in their orbital and spin indices, and are constructed to be traceless. 

^We effectively assume that in both 4-dimensions and d-dimensions, the operation of rotating tt rad 
about the z-axis is well defined, and is equivalent to flipping the direction of all axes except z-axis. 
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• Case 3: for ^P 2 , the wave functions with \Jz\ = 1 can be separated from those with 
\Jz\ = 0,2 by requiring that they are antisymmetric under flipping the direction of 
all axes except z, similar to the Case 1 above. The wave functions with = 0 and 
\Jz\ =2 can be further separated by requiring the wave function with = 0 to have 
a SO{d — 2) symmetry in the space perpendicular to 2 -axis. 

• Case 4: for ^Pi, the wave functions with \Jz\ = 1 can be singled out by requiring that 
they are antisymmetric under flipping the direction of all axes except z, similar to 
the Case 1 above. We define the symmetric part under this operation as the wave 
function with Jz = 0.^ 

Notice that in d-dimensions, there is no unique way to group states into categories with 
different J and \Jz\- Different grouping methods are equally good as long as they are 
consistent and give the correct decomposition at d = 4. They could serve as different 
schemes in the dimensional regularization. 

As an example, we consider the polarized LDMEs for J/ijj production below, which 
is of great phenomenological interest, while giving definitions of all 5-wave and P-wave 
NRQCD LDMEs in Appendix A. The most important [QQ{n)] channels in Eq. (2.15) for 
J/V' production are ^5^^, and ^pf^ (see Table. 1). For channel, the wave 

function for polarization-summed QQ-pair in its rest frame can be chosen as a (d — 1) 
dimensional vector^ 


I^QQpSi]) = (3.1) 

where ^|J (y) is the annihilation (the creation) operator for a heavy quark (a heavy antiquark) 
in the center of mass frame of the pair, with j = 1... (d — 1) are the basis 

states/vectors of the Hilbert space for and cP are the corresponding coordinates. 

In this Hilbert space with d = 4, the 3x3 matrix representation of the flipping operation 
described in the Case 1 above is diag{—l, —1,1), which can be easily generalized to (d — 
1) X (d — 1) matrix representation: 

Dfj = dm5(-1,--- ,-1,1), (3.2) 

where ellipsis represents (d —4) copies of —1. With Eq. (3.2), we separate the “longitudinal 
polarization” space spanned by the basis state {'0V^y;|O)} from the “transverse polarization” 
space spanned by the basis states {il^ with includes all values of j defined in 

^Different from the 4-dimension case, wave functions with \Jz\ = 0,1,2 for ^Pi channel are allowed in 
d-dimensions. Wave functions with \Jz\ = 1 are transversely polarized at d = 4. Wave functions with 
|J^| =2 have different parity from wave functions with \Jz\ = 1, and vanish at d = 4. It is then natural to 
consider | Jz| =2 wave functions as longitudinally polarized in d-dimensions. 

“^The coupling of two spinors in d-dimensions also have high-spin representations. On the one hand, their 
contribution to spin-0 and spin-1 representation is at least suppressed by because of the heavy quark 
spin symmetry [35]. On the other hand, one could always ignore them by using a different dimensional 
regularization scheme, where spinors are in 4-dimensions. 

Here we only consider the components at leading powers in v, relativistic corrections can be constructed 
by inserting covariant derivatives, as in Ref. [4]. 
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Eq. (3.1), except j = z. We could then construct the definitions of polarized NRQCD 
LDMEs for in d-dimensions as 

(3.3a) 

(3.3b) 

where the summation of j± is suppressed, an^ is the annihilation operator of heavy quarko- 
nium H with polarization A, and (^^_ 2 )( 2 N ) ^ normalization factor. The corresponding 

NRQCD LDMEs for channel can be obtained by removing 1/(2W) hr Eq. (3.3) and 
inserting one SU(N) color matrices T“ to the operator of each heavy quark pair. 

For J/'ip production, the analysis of the ^Pf’^ channel is greatly simplified if we assume 
the heavy quark spin symmetry, as explained at the end of Subsection 2.3. Since the 
spin of the produced J/ij} must be the same as the spin of the intermediate QQ-pair in the 
leading power of v, the orbital angular momentum of the heavy quark pair can be effectively 
summed over. In this way, we can apply the same flipping operation as in Eq. (3.2), but 
only to the spin index, and we obtain 

OH>(S,Tp[8|) ^ (3.4a) 

OH8(3,lp[8|) ^ (3.4b) 

where 

= ^HDx) - {D^)^X- ( 3 . 5 ) 

A complete treatment of the channel without summing over the orbital angular mo¬ 
mentum needs QQ channels offdiagonal in {n,h), as analyzed in Section 2. For simplicity, 
we do not consider individual offdiagonal LDMEs as their overall contributions can be 
included in the LDMEs in Eq. (3.4). 

4 NRQCD factorization on input FFs 

The NRQCD factorization formula derived in Eq. (2.15) can also be applied to FFs at the 
initial scale //q- To be specific, the input single parton FFs and the QQ-paii FFs to a 
polarized heavy quarkonium can be written in the following factorized form^ 

Df^H^,{z;rnQ, fio) = (z; mg, /tq, A'-a)(a^a))? (4.1a) 

'^\QQ{K)]^Hy^i (.2) Cl) C2) '^Qi ^ ^ Cl) C2) IHQj A^Oj ’ 

n\ 

(4.1b) 

®In this paper we work in the hadron helicity frame, in which the mixed fragmentation with a single 
parton in the amplitude and a heavy quark pair in the complex conjugate of the amplitude does not exist, 
following the same argument in Section IIB of Ref. [8] 
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where > 2mQ is the QCD factorization scale, ~ mg is the NRQCD factorization 
scale, and A (A') denotes the polarization of the intermediate NRQCD QQ pair (observed 
heavy quarkonium H). In Eq. (4.1a), / = Q,Q,q,q, g is the flavor of the fragmenting 
parton, and the variable z is the light-cone momentum fraction of the parton taken by the 
quarkonium H. In Eq. (4.1b), k = or represents the vector, axial-vector or 

tensor spin states of the fragmenting heavy quark pair, respectively, where the superscript 
labels the singlet (1) or octet (8) color state of the pair. The variables Q and ((2 in Eq. (4.1b), 
dehned in Refs. [19, 20], represent the relative light-cone momentum fractions of the heavy 
quark pair in the amplitude and its complex conjugate, respectively. 

The short-distance coefficients (SDCs), d’s in Eq. (4.1), are insensitive to any particular 
quarkonium state Hy, which is an immediate consequence of the factorization. We can 
derive these SDCs by replacing the hnal heavy quarkonium state by an asymptotic heavy 
quark pair [QQ{n'y)], 


(4.2a) 

'^[QQ]{K)^[QQ{n')]iz,Ci,C2;mQ) = diQQ](^)^[QQ(nx)]i^^ Cl, C 2 ; 


[QQ(nx)] 


(4.2b) 


n.\ 


where the dependence on the factorization scales is suppressed. In Eq. (4.2), the LHS can 
be calculated with perturbative QCD, while the LDMEs on the RHS can be calculated with 
perturbative NRQCD. By matching the LHS and the RHS, the SDCs can be obtained order 
by order in ag¬ 
in this paper, we follow the convention used in Refs. [19, 20], and expand the SDCs as 


[z] mq, go) = ^ 


n\ 


+ ( 


TT / f^lQQi'^x)] 


(z; mg, ^ 0 ,/^a) + ©(a^)} x 




m 


2L+3 

Q 


(4.3a) 


'^[QQ{k 


Cl) C2, /^o; mg) ^ { 4 QQ(K)]^.[QQ(nA)] ‘’1’ ^2, go', mQ,gA) 


n\ 




7 [QQ{f^)]^[QQ{‘^\)] 


(z, Cl, C2, go-, mq, g^) + 0 (aj ^ x 




[QQln^)] 


M) 


m 


2L+1 

Q 


(4.3b) 


Now with the d-dimensional polarized NRQCD LDMEs and the projection operators 
derived from these LDMEs, we can apply CDR to regularize all the divergences in the cal¬ 
culation of the LHS of Eq. (4.2), with a particular polarized \QQ{n'y)] state. In the NLO 
perturbative calculation of the LHS, there exists (1) UV divergences, which are removed by 
QCD renormalization and the renormalization of composite operators; (2) rapidity diver¬ 
gence, which are canceled when adding up all Feynman diagrams; (3) Coulomb divergence 
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and IR divergence, which are canceled between the LHS and the NRQCD LDMEs on the 
RHS, if NRQCD factorization is valid. The cancellation of the Coulomb divergence and IR 
divergence needs NLO calculation of polarized NRQCD LDMEs, which are also given in 
Appendix A. 

These polarized SDCs and the polarization-summed SDCs are not independent. With 
the definitions of the polarized NRQCD LDMEs in Appendix A, one can obtain the relation 


?(1 or 2) _ 1 ^ t(l or 2) 

W^[Q<9(n)] 

^ A 

j(0 or 1) _ 1 

^lQQM]^[QQ{n)] ~ Z^®[QO(K)]^[QQ(nA)]’ 
^ A 


(4.4a) 


(4.4b) 


where di^ m and m are the polarization-summed SDCs, and are 

f^lQQin)] [QQ(«;)]-s>[QQ(n)] ^ ’ "■ 

number of polarization states for the channel n, both of which can be found in Refs. [19, 20]. 
The summation of A runs over all polarization states of n. In Appendix B and C, we only 


;(0 or 1) 


list the d 


(1 or 2) 


and d 


(0 or 1) 


^f^lQQin,)] — “[QQ(^)]^[QQ(nD] polarization state for ^Pi, and ^Pi 

channels, and two polarization states for ^P 2 channel. One can derive the rest by using the 
relations in Eq. (4.4). 


5 Discussion and Summary 

In this paper, we derived all nonvanishing channels of producing a heavy-quark pair (n, h) 
for polarized heavy quarkonium production in NRQCD factorization. We introduced a 
scheme to generalize the polarized NRQCD LDMEs to arbitrary d-dimensions. With these 
d-dimensional NRQCD LDMEs, we apply CDR to calculate the single-parton and heavy- 
quark pair EEs, up to the NLO in ag, for the production of all polarized S-wave and 
P-wave heavy quarkonium at the input scale /io P 2mQ. We find that all perturbative 
divergences are canceled at this order, and all derived SDCs are finite as expected. With 
our results, as well as the evolution kernels and the hard parts calculated in Refs. [8, 18], 
the QCD factorization formalism can be used to generate predictions for the polarization 
of produced heavy quarkonia in the hadron helicity frame, which have been measured at 
the Tevatron and the LHC. 

In addition to the surprisingly large NLP contribution to the yield of J/'0 [23], our 
results in this paper further demonstrate the potential importance of the NLP to the J/'0 
polarization. In Table 2, we show the LO contributions from different channels to the J/'0 
polarization. The single-quark fragmentations are not important since they are suppressed 
at large z [19]. The contribution with a fragmenting QQ pair in a tensor state is suppressed 
by the hard part [18]. These suppressed channels are not showed in Table 2. For complete- 
ness, we also list the Wq channel, which contributes to unpolarized J/ijj production. 

Contrary to the contribution at the LP, most production channels at the NLP contribute 
to longitudinally polarized J/'ijj. Recall that the NLP contributes to the production rate, 
both directly by the NLP term in the factorized cross section, and indirectly through the 
mixed kernels in the evolution equations of the single-parton EEs (see Ref. [8] for more 
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Table 2. Contributions of LO FFs to the J/'0 polarization. The labels “T”, “L”, and “Un” represent 
transversely polarized, longitudinally polarized, and unpolarized J/'f/’j respectively. 



35 W 

3^|8] 

3pM 

l^[8] 
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T 



^[1] 

L 




til®] 


L 

L 


all] 





a[8] 



T 

Un 


details). Especially, the indirect contribution can increase the longitudinal component of 
produced J/ij) from gluon fragmentation when the gluon FF is evolved from the input scale 
Ho > 2mQ to the hard scale fi ^ Pt- Although the details need more studies, we believe 
our results in this paper will be very helpful for understanding the polarization of heavy 
quarkonium production. 
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A Polarized NRQCD LDMEs 

This appendix is organized as follows. In Subsection A.l, we list the definitions of all 
essential d-dimensional polarized NRQGD LDMEs with the 5-wave and P-wave heavy- 
quark pair. In Subsection A.2, we present the projection operators for the calculation of 
SDGs. In Subsection A.3, we expand the polarized NRQGD LDMEs to NLO in which 
are necessary for the full cancellation of IR divergences in the NLO calculation of SDGs. 

A.l NRQCD LDMEs 

With the method expained in Section 3, we give our definitions of normalized NRQGD four- 
fermion operators for polarized heavy quarkonium production in an arbitrary dimension d. 

r“x, (A.la) 

(A.lb) 
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% ^ 


t 


(d-2)" ^ 2^ V 2 


OH,(ip|sj_) ^ yii-L^‘)T‘i.(a^„a„,)i,H-]fS‘)T‘x, 


I ^ 


1 '^7^, 


(A.lc) 

(A.ld) 


O^-^pf},) = cT^l)r“V^(al,^aHjV^n-2'^ (A.le) 


1 


0j/.(3p[81 ) ^ zj 

(A.lf) 


(nH\ /Spls] N __ 

Ly ^2,tt) ~ (j _ 


{d-l){d-2) - 2 


X 


t(_l(l^ {A ^fcx} _ 
^ 2^2 


SJi 


d-2 


^T ■ crT))r“V' 


- JZT^'^T ■ crT))r“x, 


(A.lg) 




1 


X\--D ^^^a^^)T‘^^Pia^^^aH,)^PH--D (A.ih) 


d- 1 ^ ' 2^ d-2 

- jzn^ 


Dt ■ o-t)) T^X, 


Hx(3p[S]\_ 1 


x^(-^^ • 0-) r“V'(oL^aHjV’^(-^^ • 0-) r“x, 
o^M' 4®') = x^r“V’(al,^apjv^tr“x, 


/'opLoj^ _ 

^ d-1 


(A.li) 

(A.lj) 

(A.lk) 


where the summations of j± and k± run over all directions perpendicular to i-axis, and 

1-^ 


Dx = 'iPHDx) - {DipYx, 

D X • ctt ^ D • ct — D , 

^ [A ^fcx _ ^ix^ 

^ lix ^fcxl = ^ix ^kx_ ^fcx ^ix. 


(A.2a) 

(A.2b) 

(A.2c) 

(A.2d) 


Subscripts TT, T, and L represent the non-relativistic QQ pair with | J^j = 2, | J^l = 1 and 
Jz = 0, respectively. Polarization of the heavy quarkonium H is labelled by A, which can 
be TT, T and L, depending on the specihc state of the H. For completeness, we also 
list the operators for states ^Pq and ^Sq, which are unpolarized. 

As mentioned at the end of Section 3, the ^Pf^ channel can be greatly simplihed by using 
the heavy quark spin symmetry, so that we only need two NRQCD LDMEs in Eq. (3.4): 


1 


0 H,(3,Tp[8]) ^ (A.3a) 


1 


QHx^3,Lp[8]-^ _ 




I 


(A.3b) 


The color-singlet operators can be obtained from their color-octet counterparts in Eqs. (A.l) 
and (A.3) by removing the two T“’s and multiplying the factor l/(2Aic). 
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The LDMEs defined in Eqs. (A.l) and (A.3) have been normalized by the number of 
spin states, 



= N^pI^ = = d - 2, 

(A.4a) 


- AiPi,^ - Asp^ ^ - 1, 

(A.4b) 

^^Pt 

= ^(d-2)(d-3). 

(A.4c) 

A/NR 

Sp^TT 

= l(d-l)(d-2 )-l. 

(A.4d) 

;VNR 

^^3,Tp 

= id-l)id-2), 

(A.4e) 

-‘''s.np 

= d- 1. 

(A.4f) 


By adding the number of spin states with the same J but different \Jz\-, we can retrieve the 
normalization factor for unpolarized heavy-quark-pair in Ref. [19, 20]. With the definitions 
above, it is straightforward to check that adding up all operators with the same J but dif¬ 
ferent \Jz\ weighted by the number of spin states, we reproduce the conventional definitions 
of unpolarized NRQCD LDMEs at arbitrary d dimension, 

(A.5) 

ny 

which correspond to definitions in Ref. [4] by setting d = 4. ® 

A.2 Projection Operators 

With the definitions of d-dimensional polarized NRQCD LDMEs, we derive the projection 


operators in the same way as in the polarization-summed case, 

P^sl,=P^pl,=^Tip), (A.6a) 

P^sl^=Prpl^=Ffip), (A.6b) 

Pij.1^ = ^ [nP{p)FfP{p)+Ff{p)Ff\p)-Ffip)Fl^\p^ 

, (A.6c) 

^ (p)Pf' (p) - Pf' (p)Pf' ip)) , (A.6d) 

^ (Fr\p)Ff(p) + pf'(p)P“'^(p)) - ^Ff(p}F^P(p), (A.6e) 

= I (n“'(p)pf'(p) 

+pf'(p)Pj^'(p)) , (A.6f) 

^^P^L = “ j^K^(p))(P^'^'(p) - J^K'^'(p))’ (A-6g) 


®Note that our definitions of color singlet NRQCD LDMEs differ from the definitions in Ref. [4] by a 
factor of 1/(2W), while the definitions of color octet NRQCD LDMEs are the same. 
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= P““'(p)pf'(p), 

p^R ^ 


where the superscript “NR” refers to NRQCD, and 


TroaaV N aa' , P ^ 

™± [P) = —g H-:-^- TTi^n n 


p ■ n 


{p ■ h) 


^aa',^ PP p“n“+p“n“, P^ 

\P) ~ —9 -^-^ ^- TTi^n n 


p ■ n 


{p ■ n) 


Q. d 

P““'(p) = Pff“'(p) + Pi“'(p) = -5"“' + ^ ^ 


2 ■ 


(A.6h) 

(A.6i) 

(A.6j) 

(A.6k) 


(A. 7a) 
(A.7b) 


(A.7c) 


a and j3 {a' and fj') are the indices for the orbital angular momentum and spin of the 
heavy quark pair in the amplitude (complex conjugate of the amplitude), respectively. By 
adding up all projection operators with different | J^| but the same J, we can get the same 
unpolarized projection operators in Refs. [19, 20]. By setting d = 4 for the with 

A = L,T,TT, we retrieve the results in Ref. [37]. 


A.3 Expand the polarized LDMEs to NLO in perturbative NRQCD 

To complete the cancelation of all IR divergences in the NLO calculation of SDCs, we need 
to expand the polarized NRQCD LDMEs to NLO in powers of a^, the same as what was 
done in the polarization-summed case. We calculate the NLO NRQCD corrections for four- 
fermion operators using similar method as described in Ref. [38]. For our purpose, we only 
give the NLO correction of 4-fermion S-wave operators: 

- Ce^\^{d - 2)(C)[W]A(lp[8^)^LO ^ 

^ (C>[WA(l5[8lj)LO 

- C,[Cf ((d - 2)(0['3Q]A(ip[i^))LO ^ ^0[Q(9]A(ip[il.))LO^ 

+ Bf ((d - 2)(C>[WA(lpj8l^))LO ^ ^0[OQ]A(lpj8l.)^LO^ 
(0 [QQ]a(3^[i^)^ ^ wP, (A.8c) 

(0[QQ]a(35[1^)^ ^ (0[WA(35[11.)^L0 _ 

^0 [QQ]a(3 ^[8^)) = (0 [W]a(3 ^[8 ]^^^lo _ c.I^Cf wP + Bf LfJ?'}, (A.8 e) 

^0[QQ]a(3^[8|^)^ ^ ^0[QO]a(3^[8]. )^LO _ c,^Cf + Bf Lf]®'}, (A.8 f) 
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where Cp = Bp = 


and C, = TT^7^(47re-T'^)" I ^ 


and are defined 


STrmg e/fl 


MA 


In Eq. (A.8), 


as 


wP ={d - i)(^o^QQh(3,Tplb]^^LO _^ ... 




d- 


d-3 


+ li^(o[«oi.(3p[i.| ))L0 + (l;:i)^L^(oie«i3(3pi‘i,^)) 

1 


3pM 


+ + (o[QQ]x(3pib]^^'^LO ^ . . . 


(d- l)(d-2) 


(A.9a) 


wP ={d - i)(C)[QQ]A(3,ip[b]))Lo . 


1 


d- 1 
d- 2 

+ 


(^q[QQ]\(^ 3p^6] ^ ^ LO 


+ 


d- 2 


0 - '2 
3pW \\LO . d -2 


(^Q[QQU^3p[b]^^^LO 


1,TJ 


(A.9b) 


^0[QO]A(3pjfe^^^LO p - tp[QQ]x(3p[b]^^^'^LO 


+ 


where {O^QQ]^{^Ppib]))^o ^^[QQIa(3 ,Lp[fe])^LO are the leading order terms of LDMEs of 
NRQCD operators in Eqs (A.3), respectively, and ellipsis denotes terms that are irrelevant 
for our calculation. 

To derive the expressions in Eqs. (A.8), we replace d by (4 — 2e) only for the dimension 
of loop momentum in the calculation while keeping d’s in other places untouched, and 
define a renormalization scheme for LDMEs to subtract all terms that are proportional to 

, which is the same as by replacing ejp by ejjv Our renormalization scheme is 
similar to the conventional MS scheme, but not exactly the same. The advantage of our 
renormalization scheme is that, if LDMEs are dehned differently by a factor of 1 + 0(e), 
the d dependence in Eqs. (A.8) will be modihed correspondingly, but the calculated hnite 
short-distance coefficients given in the next two Appendixes are unaltered . 


B Single-Parton Fragmentation Functions to a Polarized Heavy Quarko- 
nium 


In this appendix we list the SDCs for all single-parton fragmentation functions to S-wave 
and P-wave polarized QQ-pair up to order 0{a‘^) dehned in Eq. (4.3a). The polarized 
EEs and unpolarized EEs are related by Eq. (4.4). Therefore for outgoing QQ with n -|- 1 
polarizations, we only give n polarized EEs below. The other one can be derived by using 
Eq. (4.4) and the unpolarized EEs calculated in Ref. [19]. All channels that vanish for all 
polarizations are not listed. 


B.l Gluon FFs 
Leading order 

- f8i = 0. 

5^[QQ(3s1%)] 
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Next-to-leading order 
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^ X d^^^ 
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^ X 

Cf g^[QQ?pfl,, t)\' 


B.2 Different quark FFs 

J(2) ^ 2(1-2)^ 

q^[<3Q(3s(®^)] Ncz{ri 22 - 42 -F 4) ’ 

where 77 = nij^/mQ. 

B.3 Same quark FFs 
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(B.4) 

(B.5) 

(B. 6 ) 
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g^[Q<3(3p|)[)] 

- [11 

Q^[Q<3(3p))^)] 


-^5922^ - 6402-^320), 


{N^ - 1)^ ^(1 

15N3 (2 


2f 


(352® 


3122® -f 11362^ 


- 2 OI 62 ® -h 18722^ - 9602 -h 320), 


(AT,^ - 1)^ 2(1 - 2 )^ e 

15N3 ( 2 - 2 ) 8 ^ 


5802® -h 16282"^ 


(B.IO) 


(B.ll) 


(B.12) 
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- Ill 

Q^[Q<9(3S'%)] 


Q^[QO(ip|®[)] 

Q^[Q0(3p|®[)] 


d 


( 2 ) 

''Q-^[QQ(3P 


[8] 


2,T or TT) 


- 1872^3 + 1328^2 - 512z + 128), 

(5z^ - 32z^ + 68z^ 


8{N^-1)^ zil-z)\^, oo-,3 , «o-,2 


15A^| {z - 2)8 
67V| i(z - 2)6 " 12iVc22(^ 


+ z2(3z^ - 24z^ + 64^2 - 32z + 16) , 
_ \ 

(A^|_1)2 Q^[QQ(lp[)[)]’ 

_ \ _J(2) 

(7V2_1)2 Q^[QQ(3p[)[)]’ 

_i_J(2) 

(A^l _ 1)2 Q^[QQ(3pJ8^ 



(B.13) 

-32Z + 16), 

(B.14) 

1 

1 

to 

to 

(B.15) 


(B.16) 


(B.17) 


(B.18) 


B.4 8p^ operators with orbital angular momentum summed 

Here, we list the SDCs corresponding to the NRQCD LDMEs o'" §] dehned in Eq. (A. 3). 


d 


( 2 ) 

g^[(3Q(3Pp[l])] 


9^[<9Q(3’-^pM)] 

Q^[QQ{^’^PW)] 


d 


( 2 ) 

Q^lQQCS.i'PM)] 


1 


2{z - l)(z^ + 8z - 12) ln(l - z) 


2NcZ^ 

+ zi2z^ + z‘^ - 28z + 24) 

Bp ?(2) 

Cp ^ 9^[Q0(3'^pm)] 

~ y ~ (23z^ - 1922*^ + 676z^ - 1120^^ 

6iV3 (_2 _ 2)8 ^ 

+ 1104^^-512^ + 192), 

1_h(2) 


(^■2 _ l)2“'Q^[Q0(3’-^Pm)]' 


(B.19) 

(B.20) 


(B.21) 

(B.22) 


B.5 Comparing with other calculations in the literature 

Polarized EEs in g —)• QQ(^S^^j^) channel has been calculated in Refs. [39, 40], our result in 
Eq. (B.5) conhrms their results. 

EEs from g —)• QQ{^P\^^) and g —)• QQ{^P^^) channels in different polarization states 
have been calculated in Ref. [41] with cutoff regularization scheme. Comparing their re¬ 
sults with ours in Eqs. (B.2), (B.3) and (B.4), we hnd that if we change log{g,\/mQ) by 
\og{g.\/mq) + 5/3 in our results, and take into account the different normalization of the 
LDMEs, we can reproduce their results in Ref. [41]. The 5/3 difference is due to different 
regularization scheme and was also realized in the calculation of polarization-summed EEs 
by authors of Ref. [42]. Therefore, our results of EEs from g —)■ (5(5(^P|^^) and g —)■ QQ{^P^^ ) 
channels in different polarization states are consistent with the results in Ref. [41]. 

Polarized EEs in g —)■ QQ(^5'pj;^), Q —)■ (5(5(^5'[^{^) and Q —)■ (50(^‘S'fp) channels have 
been reported in Ref. [36] when we are preparing our paper. ^ The authors of Ref. [36] 
present their results with NRQCD LDMEs dehned slightly different from ours. After taking 

^Our results were published previously in Ref. [43], in which there was a typo for q —>■ (5(3(^S'[*/) channel. 
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into account the difference of the LDMEs, our results in Eqs. (B.9), (B.IO) and (B.15) are 
consistent with theirs. 


C Double-Parton Fragmentation Functions to a Polarized Heavy Quarko- 
nium 

In this appendix we list the SDCs, up to order ©(a^), for all heavy quark pair fragmentation 
functions to a polarized QQ-pair in S-wave and P-wave, as defined in Eq. (4.3b). The 
polarized EEs and unpolarized EEs are related by Eq. (4.4). Therefore for outgoing QQ 
with n + 1 polarizations, we only give n polarized EEs below. The other one can be derived 
by using Eq. (4.4) and the unpolarized EEs calculated in Refs. [19, 20]. All channels that 
vanish all polarizations are not listed. 

C.l A-functions 

In our results below, we use the same definitions of A-functions as that introduced in 
Ref. [20]. We repeat these definitions below for readers’ convenience and the completeness. 


Ao = 4 5(Ci) 5(C2), (C.l) 

A^^ = Az^6'{Ci)S\C2), (C.2) 

= 4 [5(1 -z + Ci)± 6 {l-z- Cl)] [5(1 -z + C 2 )± 6 { 1 -z- ( 2 )] , (C.3) 

a[^' = -iz\ [5'(1 - z + Cl) ± 5'(1 - 2 ; - Cl)] [5(1 - z + C 2 ) ± S{1 - z - C 2 )] 

^ (C.4) 

+ [5(l-z + Ci)±5(l-z-Ci)] [5'(1 -z + C2)±5'(1-^-C2)] }, 

A^f = 4z^ + [s'{i-z + C2)±5'{l-z-C2)\ , (C.5) 

Alt®' = 4| {Nl - 2) [5(1 - z + Ci)5(l -z + C 2 ) + 5{1-z- Ci) 5(1 - z - C 2 )] 

^ (C.6) 

T 2 [5(1 - z + Ci)5(l - z - C 2 ) + 5(1 - z - Ci)5(l -z + C2)]\, 


Alt®’' = -4 z[{Nl - 2) [5'(1 - z + Ci)5(l - z + C 2 ) + 5(1 - z + Ci)5'(l - z + C 2 ) 

+ 5'(1 - z - Ci)5(l - ^ - C 2 ) + 5(1 - z - Ci)5'(l - z - C 2 )] 

T 2[5'(1 - z + Ci)5(l - z - C 2 ) + 5(1 - z + Ci)5'(l - z - C 2 ) 

+ 5^(1 — z — Ci)5(l — 2 ; + C 2 ) + 5(1 — z — Ci)5'(l — z + C 2 )] |) 

Alt®’" = 4^2 I(iv2 _ 2 ) [5'(i _ ^ + Ci)5'(l - z + C 2 ) + 5'(1 - z - Ci)5'(l - z - C 2 )] 

^ (C.8) 

T 2 [5'(1 - z + Ci)5'(l - z - C 2 ) + 5'(1 - z - Ci)5'(l - z + C 2 )] }, 

All these A-functions are invariant under the transformation (Ci —t —Ci) C 2 —C 2 ) and 

the exchange Cl t-7> C 2 ) including the crossing exchange (Ci —t —( 2 , C 2 —Cl)- In addition, 
Aq, a’^’, a’I];’ and a’^’ are even in both Ci and C 2 ) while Aq, a’]';’, a’]';’ and a’^’ are odd 
in both Cl and C 2 - Under the integration of Cl and C 2 with a well behaved test function, the 
asymptotic behaviors of these A-functions at z —>■ 1 are 

lim A’t’ = 0[1], lim a’^’ = 0[(1 — 2 :)^], 

z—s^l z—^1 
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Therefore, 


lim a[J^ = 0[(1 — z)], 

lim A 


2 :—>-1 

lim a[J^ =0[(1 — 2 )^], 

lim A 


2^1 

lim a!?^ = 0[1], 

lim A^ 

z—>-l 

lim a!?^ = 0[1], 

z—>-l 

Z^l 


= 0[(l-z)], 

= 0 [ 1 ], 

= 0[(l-z)], 


aW aW' aW" A[f]' 

(1 -z)’ (1-z)’ (1-z)’ (1 -z) 

do not exhibit any pole at 2 : = 1. 

C.2 Leading Order 


[Qg(^M)]^[Q(3(3p'i4)] 

[QQ(a[il)]^[QQ(ip|^^)] 

[(3(3(a[il)]^[QQ(3p[i^)] 

d^^'> 

d^°^ 

[QQitm^[QQ(^Pi%T)] 

[QQ(tM)]^[QQ(3p[i[)] 


0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 . 


The corresponding color octet channels also vanish. 


C.3 NLO - Vector 


r, - Ml =0, 

d^^\ 

[QQivW)]^[QQePi^l)] 




[QQivm^lQQi^Pi'lT)] 


= 0 , 
= 0 , 


d 


( 1 ) 

[qO(^I'')]^[<30(3s1%)] 8N, 


2,TTJ 

[ 8 ] 


^ -{l-z)zA^^\ 


(C.9) 

(C.IO) 

(C.ll) 

(C.12) 

(C.13) 

(C.14) 

(C.15) 

(C.16) 

(C.17) 

(C.18) 


(C.19) 

(C.20) 

(C.21) 

(C.22) 
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( 1 ) 


1 


[QQ(i,M)]^[QQ(ip|%)] 12N, 

( 1 ) ^ 
[QQ(i,M)]^[QQ(3p[«^)] gQNc 

.W 


{l-z)zA^l\ 

|8Ao(5(l-z)(- ln[^] +21n2 + ^) 


(C.23) 


"Q 


+ 6Al^J z{l-z)-3A^f z{2z-3) 


d 


+ 4:ZA 


1 


[ 1 ] 


- 2z + 3 


}. 


( 1 ) _ _ 

[QQ(i,M)]^[QQ(3pP^)] 96iV, 


(1 - z) + 

2 

|8Ao(5(l — z)(^ — ln[-^] + 21n2 + 


(C.24) 


Q 


+ 6AW'' z{l-z)- 3AW' z {2z - 1) 


+ AzA 


[ 1 ] 




[QQ(^[ll)]^[QQ(3p[8^^)] 24A( 


(1 - z)+ 

■z(1-z)aW 


-2z + l 


}■ 


+ > 


2,TTJ 

ylS] 


1 


( 1 ) 

[Q0(^^t»])]^[Q0(3sl%)] 8Nc{N^ - 1) 

( 1 ) _ 1 


{1-z)zA^3, 


{l-z)zA^^\ 

2 

[Q0(^'f«])]^[Q0(3pl%)] “ 96Ac(A2 - 1) " 4 )Ao5(1 - z)( - ln[^ 

+ 21n 2 + + 6A^^'' z{l-z)- 3A^’' z {2z - 3) 


[QQ(i,[31)]^[Q0(iPl%)] 12Ac(A2 - 1) 

( 1 ) _ 1 


(0.25) 

(0.26) 

(0.27) 

(0.28) 


+ 42;A 


[ 8 ] 


- 22; + 3 




HI) 


(1 - z)+ 

{4(a2-4)Ao5(1-z)(- ln[ 


(0.29) 


I^A 


[QQ(i,[31)]^[QQ(3p[«^)] 96 Ac(A2 - 1) 

+ 21n 2 + + 6A^^" z{l-z)- 3A^^' z {2z - 1) 


Q 


+ AzA 


[ 8 ] 


-2z + l 




[Q0(^[®>)]^[Q0(3prW)l 24Ac(A2 - 1) 


(1 - z)+ 

z{l-z)A^l\ 


}. 


d 


( 1 ) 


= 


[gQ(,;[81)]^[gg(2S+iLMj] [Qg(,;Iil)]^[QQ(2S+ipMj]’ 

where A is the polarization of the outgoing heavy quark pair. 

C.4 NLO - Axial-vector 


(0.30) 

(0.31) 

(0.32) 


M, - , [11 = 0> 

[Qg(aIil)]^[Q(3(iPiW )] 

d^^^- ,, - rii =0, 

[Qg(aIil)]^[Q(3(3PjW )] 

d^^\ -18] = — z)A[,*;^, 


(0.33) 

(0.34) 

(0.35) 
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?( 1 ) 


[QQ(a>il)]^[Q<3(ipl®i,)] 2AN, 


+ 2zA 


J^{4A„{(1-^)(-1„[£|]+21„2+1) 

}. 


3 

— -5; 


d 


( 1 ) ^_ 

[QQ(aIil)]^[Q(5(3p[%)] 96A, 


(1-Z)+ 2 J 

^z(l - z)(6AW'' + 9AW' + 16aW), 


d 


( 1 ) 


[QQ(aIil)]^[QQ(3pJ®^)] 


^z(l - z)(6AW" + 3AW' + 8AW), 


d 


1 


( 1 ) __ 

[QQ(a[il)]^[(3Q(3pfl )] 24A, 


( 1 ) 

[Q0(a>®])]^[Q<3(3s'%)] 

l(l)_ 

[QQ(a[31)]^[QQ(ip'®i,)] 


■z(l -z)aW, 


d 


1 


8Ae(A2 - 1) 
1 


z(l-z)A'f], 


24A,(A2-1){^^^" 4)Ao(5(1 z)( 

' }. 


Q 


+ 2 In 2 + —^ -j- 2 z A[j_^ 


3 

— -5; 


L(l-z)+ 2 J 


d 


[Q0(af81)]^[Q<3(3pl®i,)] 

( 1 ) 

[Q0(aM)]^[QQ(3pi%)] 


1 


d 


( 1 ) 


96Ac(A2 - 1) 

1 

96Ac(A2 - 1) 

1 


[Q0(a[il)]^[QQ(3p>*4)] 24 Ac(A2 - 1)' 

= d^^\ 

lQQ(alS])]^[QQ(2S+lLlj[)] [Q0(aIll)]^[QQ(2S+iLMj]’ 


z{l - + 9A^'' + 16A^'), 

z{l - z) (GA^!''' + 3A^'' + 8A^^), 
z{l-z)A^^\ 




(C.36) 

(C.37) 

(C.38) 

(C.39) 

(C.40) 

(C.41) 

(C.42) 

(C.43) 

(C.44) 

(C.45) 


where A is the polarization of the outgoing heavy quark pair. 


C.5 NLO - Tensor 




= 0 , 


r, - Ml =0, 


d 


( 1 ) 


[QQ(tM)]^[QQ(3p[i[)] 


= 0 , 


= 0 , 


[Q(3(hil)]^[Q(3(3pli4)] 

d^^^- r, - Ml =0, 
[QQ(tM)]^[QQ(3pW )] 


( 1 ) ^ _ 

[Q(5(tM)]^[QQ(35f^)] 16 A, 


—z(1-z)aW, 


(1) ^ _L 

[QQ(tM)]^[Qg(ip'«[)] 48Ac 
( 1 ) _ 1 


[QQ(tM)]^[QQ(3pf[)] 48Ac 


z(1-z)aW, 
2Ao5(l - z) 



(C.46) 

(C.47) 

(C.48) 

(C.49) 

(C.50) 

(C.51) 

(C.52) 
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+ 2 A 


[ 1 ] 


-42 + 1 


}■ 


d<'l 


241V, 


( 1 - 2 ). 

+ {2A„4(l-^)(-ln[£|-]+21„2 + 


+ 2A 


[ 1 ] 


-z + 1 


h 


?( 1 ) 


[Q(5(tW)]^[QQ(3p[8[)] 288N, 


( 1 - 2 ) 

.{4A„i(l-^)(-l.,[A]+21„2+l 


(C.53) 


(C.54) 


+ 12 2 AW" (1 - z) -3zA^l^' (32 - 2) 


[ 1 ]' 


+ 2 0 A 


[ 1 ] 




1 


'[QQ(t[«l)]^[QQ(3s[%)] 16iVc(A2 - 1) 

_ = __ ^ . 

[QQ(t[ 81 )]^[QQ(lp|«[)] 48 Ac(A 2 - 1 ) 

1 


. 72 + 5 

}, (C.55) 

- 2 )+ 

-02(1 - 2)A?1 

, (C.56) 

-02(1 - 2)AL®1 

, (C.57) 


[QQ(t[8])]^[QQ(3p[8[)] 481 Vc(A2 - 1)1*'^'" z)(^ ln[ 

-42 + i]}, 


i4- 

^Q' 


+ 21n2 + ^) +zAf 


?(i) 


1 


[Q<3(i[«l)]^[QQ(3^’i®W)] 24A^c(iVc^ - 1) 


?(i) 


+ 2 In 2 + —^ + 2 A^_^ 


1 


[QQ(t[8])]^[QQ(3p[8[)] 288Nc{N^ - 1) 


( 1 - 2 ) 

{(+-4)A„4(l-^)(-ln[A] 

^ ^ '''-Q 

7i-^- Z + l 

L(1-2)+ JJ 

{2(+-4)A„4(l-2)(-ln[A] 

^ ^ '''-Q 


(C.58) 


(C.59) 


+ 2In 2 + ^) + 12 2 A[f'" (1 - 2 ) - 3 2 A^^' (32 - 2) 


+ 22 A 


[ 8 ] 


1 


- 72 + 5 


}■ 


‘ ■ r [3] M ^rr4^2*rilM V r [8] m’ 


(C.60) 

(C.61) 


-(1 - z). 

'[QQ(48])]^[QQ(2S+ilWj] “ “[QQ(tW)]^[Q0(2S+iLMj]’ 
where A is the polarization of the ontgoing heavy qnark pair. 

C.6 ^Pj operators with orbital angular momentum summed 

In this snbsection, we list all SDCs corresponding to the NRQCD LDMEs 8Ap[i or 8] dehned 
in Eq. (A. 3). 

C.6.1 Leading Order 


- 0 

“[(3Q+M)]^[QQ(3.T’p[il)] 

_ — Q 


(C.62) 

(C.63) 
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— 0 


(C.64) 


The corresponding color octet channels also vanish. 

C.6.2 Next-to-leading order 


“[Q(3(t;M)]^[QQ(3>rp[i])] 

“[QQ(,;[il)]^[QQ(3,rp[8])] 


d^^'> 

“[Q(3(t;M)]^[QQ(3>rp[i])] 

d^^'> 

“[QQ(,;M)]^[QQ(3,Tp[8])] 
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( 1 ) 


“[QQ(aI81)]^[QQ(3,i^p[il)] 

“[QQ(al81)]^[QQ(3Pp[81)] 


“[QQ(hil)]^[QQ(3PpM)] 

[QQ(d8l)]^[QQ(3Pp[il)] 

“[QQ(h81)]^[QQ(3,rp[81)] 


0 , 

0 , 

0 , 


(C.65) 

(C.66) 

(C.67) 


8N, ' 

+ Af+ Af}, 


{a?1i„[£|]-a[;i 

21n(2(l -z))-l 

L TTIq 



d, 


( 1 ) 

[QQ(^[l])]^[QQ(3,Tp[8])p 
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8Nc{N^ - 1) 


Kl-^){A^n[£^] 


- A 


16A, 


[ 8 ] 


21n(2(l -z))-l 


+ A 


[ 8 ] 




z(l-z)] (aW + A^'+ AW") ln[^] 


[ 1 ]' 


.[ 1 ]' 




m, 


-21n(2(l-z))-3 


-A 


[ 1 ] 




Q 


d, 


( 1 ) 

[QQ(aW)]^lQQ(3,LplS])y 


(C.68) 

(C.69) 


(C.70) 


(C.71) 

(C.72) 


16Ac(A2 - 1) 

- 21n(2(l -z)) -3 




[8]+A[«l' + Al«]")[ln[4] 
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-A 


[ 8 ] 






d, 


( 1 ) 

[QQ(hll)]^[Q(9(3Pp[81)]> 

1 


IGNciN^ - 1) 


z(l-z)(2A[f' +Aff]' + A'f'"). 


(C.73) 

(C.74) 

(C.75) 

(C.76) 
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